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\  Abstract 

G/l  ^  X 

Under  appropriate  conditions,  ocean  waves  may  be  modeled  by 
certain  nonlinear  evolution  equations  that  admit  soli  ton  solutions  and 
can  be  solved  exactly  by  the  inverse  scattering  transform  (1ST).  The 
theory  of  these  special  equations  is  developed  in  five  lectures.  As 
physical  model s^these  equations  typically  govern  the  evolution  of 
narrow-band  packets  of  small  amplitude  waves  on  a  long  (post-linear) 
time  scale.  This  is  demonstrated  in  Lecture  I,  using  the  Korteweg- 
deVries  equation  as  an  example.  Lectures  II  and  III  develop  the  theory 
of  1ST  on  the  infinite  interval.  The  close  connection  of  aspects  of 
this  theory  to  Fourier  analysis,  to  canonical  transformations  of 
Hamiltonian  systems,  and  to  the  theory  of  analytic  functions  is 
established.  Typical  solutions,  including  soli  tons  and  radiation,  are 
discussed  as  well.  With  periodic  boundary  conditions,  the  Korteweg- 
deVries  equation  exhibits  recurrence,  as  discussed  in  Lecture  IV.  The 
fifth  lecture  emphasizes  the  deep  connection  between  evolution 
equations  solvable  by  1ST  and  PalnlevS  transcendents,  with  an 
application  to  the  Lorenz  model. 
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Part  I 

The  Physical  Meaning  of  Equations  with  Soli  tons 


One  of  the  important  recent  advances  in  mathematical  physics  has 
been  the  discovery  that  certain  nonlinear  evolution  equations  can  be 
solved  exactly  as  initial  value  problems,  using  a  method  that  may  be 
called  the  inverse  scattering  transform  (1ST).  One  of  the  remarkable 
features  of  many  of  these  special  equations  is  that  they  admit  as 
solutions  extremely  stable  objects  called  solitons.  For  convenience, 
we  may  use  "soliton  theories"  loosely  to  describe  all  of  these 
equations  solvable  by  1ST. 

It  may  be  regarded  as  something  of  a  miracle  that  there  are 
nonlinear  evolution  equations  that  are  completely  integrable.  It  is 
a  second  miracle  that  several  of  these  equations  arise  naturally  as 
models  of  various  physical  systems,  including  aspects  of  ocean  waves. 
My  lectures  at  this  School  will  describe  some  of  the  mathematical 
theory  that  has  been  developed  in  the  last  fifteen  years  to  solve 
these  special  equations.  Before  discussing  the  theory,  however,  it 
may  be  useful  to  give  some  idea  of  the  sense  in  which  these  soliton 
theories  model  ocean  waves.  Consequently,  in  this  first  lecture,  I 
will  try  to  describe  the  physical  meaning  of  equations  with  solitons. 

The  spectrum  of  ocean  waves  is  large  and  diverse,  and  soliton 
theories  describe  a  rather  small  part  of  it.  In  order  to  see  the 
context  in  which  soliton  theories  arise,  we  may  attempt  a  crude 
classification  of  ocean  waves,  admitting  in  advance  that  the 
classification  probably  will  be  incomplete. 

Let  us  first  classify  ocean  waves  on  the  basis  of  whether  the 
wave  amplitudes  are  large  or  small.  Large  amplitude  waves  may  break, 
among  other  things.  Fully  nonlinear  theories  are  needed  to  describe 
them,  and  will  be  discussed  by  other  speakers  at  this  School.  I  will 
restrict  my  attention  to  small  amplitude  waves,  because  soliton 
theories  arise  in  the  context  of  small  amplitude  waves.  These  waves 


are  not  necessarily  Infinitesimal,  but  they  are  small  enough  that  we 
may  develop  weakly  nonlinear  theories  to  describe  them.  All  of  the 
soliton  theories  of  water  waves  that  have  been  developed  to  date  have 
been  weakly  nonlinear  theories. 

Within  the  context  of  small  amplitude  waves,  we  may  classify  the 
theories  that  have  been  developed  as  either  stochastic  or  deterministic. 
Note  that  this  is  a  classification  of  the  theories  (i.e.,  mathematics) 
rather  than  the  waves  (i.e.,  physics).  One  of  the  interesting  open 
questions  in  this  subject  is  to  identify  the  physical  criteria  that 
determine  whether  or  not  a  given  system  of  small  amplitude  waves  may 
exhibit  chaotic  behavior.  I  will  return  to  this  point  in  my  fifth 
lecture.  In  the  meantime,  let  us  restrict  ourselves  to  strictly 
deterministic  theories. 

By  the  time  we've  restricted  ourselves  to  deterministic  theories 
of  small  amplitude  waves,  it  begins  to  sound  as  if  we'll  end  up  with 
the  linear  theory  of  infinitesimal  waves,  described  in  detail  by 
Lamb  (1932),  Stoker  (1957),  and  Wehausen  and  Laitone  (1960).  That 
turns  out  to  be  almost  right.  Soliton  theories  arise  as  (singular) 
perturbations  of  linear  wave  systems,  and  there  are  relations  between 
the  linear  theory  and  the  soliton  theories  in  terms  of  their  physical 
derivation,  in  terms  of  their  methods  of  solution,  and  in  terms  of 
their  solutions  themselves.  Soliton  theories  are  nonlinear,  but  its 
useful  to  keep  re-iterating  the  question,  "How  does  this  relate  to 
linear  theory?" 

The  physical  relation  between  linear  theory  and  soliton  theories 
is  this.  If  the  waves  are  infinitesimal,  the  linear  theory  gives  a 
complete  description  of  their  evolution.  If  the  waves  have  small  but 
finite  amplitude,  then  the  linear  theory  breaks  down  after  a  finite 
time,  and  nonlinear  corrections  are  needed  to  extend  the  range  of 
validity  of  the  theory  to  a  longer  time-scale.  Typically,  soliton 
theories  provide  the  nonlinear  corrections  to  render  the  linear  theory 
valid  on  a  longer  time-scale.  There  is  a  short  time-scale  on  which 
the  linear  theory  applies,  followed  by  a  longer  time-scale  on  which 
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the  soliton  theory  applies,  perhaps  followed  by  an  even  longer  time- 
scale  on  which  something  else  applies. 

What  are  these  soli  ton  theories?  Here  is  a  list  of  some  of  the 
equations  that  can  be  solved  exactly  by  1ST,  and  also  model  ocean 
waves. 

1 .  Small  amplitude  waves  propagating  in  only  one  spatial  dimension. 

(a)  The  \orteweg-deVries  (KdV)  equation, 

ut  +  6uux  +  uxxx  "  0  •  (1>1) 

governs  long  surface  (or  internal)  waves.  (Here  subscripts 
denote  partial  derivatives.)  See  Korteweg  and  deVries 
(1895),  and  Hammack  and  Segur  (1974,  1978). 

(b)  Under  other  circumstances,  the  evolution  of  "long" 
internal  waves  may  be  governed  instead  by  the  modified 
Korteweg-deVries  (mKdV)  equation 

ut  -  uJux  +  uxxx  *  0  *  (1-2) 

an  equation  due  to  8enjamin  (1967)  and  Ono  (1975), 

* h  ■  »•  i'-3’ 

or  by  other  models. 

(c)  Nearly  monochromatic  (l.e.,  narrow  band)  surface  waves 
are  governed  by  the  nonlinear  Schr'ddinger  equation, 

1«t  +  4»xx  ♦  a | 2  <|i  ■  0  ,  o  -  *1  ;  (1-4) 

(Zakharov,  1968;  Hasimoto  and  Ono,  1972;  Vuen  and  Lake, 

1975). 
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(d)  The  sine-Gordon  equation, 

$xt  *  0  sin't>  •  0-5) 

governs  the  waves  that  are  nearly  neutrally  stable  in  a 
baroclinically  unstable  system  on  a  beta-plane  earth  (Gibbon, 
James  and  Moroz,  1979). 

2.  Small  amplitude  waves  in  more  dimensions. 

(a)  The  equation  of  Kadomtsev  and  Petviashvili  (1970), 

(ut  ♦  uux  +  ouxxx)x  +  uyy  =  0  ,  (1.6) 

is  a  two-dimensional  generalization  of  the  KdV  equation. 

Another  general ization  of  KdV  (in  a  different  limit)  is  the 
cylindrical  KdV  equation, 

qt  +  (2t )*1  q  -  6qqx  +  qxxx  «  0  ,  (1.7) 

(b)  The  resonant  interaction  of  three  (narrow-band)  packets  of 
nearly  resonant  internal  waves  satisfy 

Oiai  *  i y »  a*  a|  , 

0j  82  •  iy2  a*  a*  ,  (1.8) 

Dj  as  *  if*  a*  a*  , 

where  Dj  *  3t  +  (Cj  •  V)  . 

Clearly  there  is  no  point  in  trying  to  derive  all  of  these 
equations.  Rather,  let  us  derive  the  KdV  equation  fairly  carefully  by 
a  multiple  time-scale  argument,  and  I  simply  will  assert  that  the  other 
equations  may  be  derived  by  multiple  time-scale  arguments  as  well.  Me 
begin  with  the  following  "exact"  problem.  Me  assume  that  the  fluid  Is 
homogeneous,  incompressible  and  Inriscid.  It  is  subject  to  a  constant. 
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vertical  gravitational  force  (g),  and  rests  on  a  horizontal  impermeable 

bed  at  z  *  -  h  (see  Figure  1).  The  pressure  vanishes  at  the  free  Fg. 

surface,  defined  by  z  *  n  (x,  y,  t),  where  surface  tension  may  be 

acting  or  not.  (In  this  derivation  we  will  omit  it,  but  almost  no 

qualitative  changes  are  needed  if  it  is  included.)  The  motion  of  the 

fluid  under  these  forces  is  assumed  to  be  irrotational,  two-dimensional 

(3y  =  0),  and  either  to  vanish  as  1  x |  ■*  »  or  to  be  periodic  in  x 

for  all  time. 

The  governing  equations  under  these  conditions  are  well  known 
(e.g..  Stoker,  1957): 


v2q> 

=  0 

-h  <  2  <  n(x,t) 

(1.9a) 

W  = 

3<t>/3z  »  0 

2  *  -  h 

(1.9b) 

\  + 

unx  ■  w 

i 

(1.9c) 

♦t  + 

»a  1  V«>  |2  +  gn  =  constant 

j  Z  *  n(x,t)  , 

(1.9d) 

and  either 

|v<t>|  ,  n  ♦  0  [ x |  ♦  ®  (1 ,9e) 

or 

I vq> |  ,  n  periodic  in  x  with  period  L  . 

These  equations,  along  with  appropriate  initial  conditions,  uniquely 
determine  the  fluid  motion,  at  least  for  some  finite  time.  A 
consequence  of  these  equations  is  that  there  are  three  globally 
conserved  quantities: 
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M  = 


net  mass  of  wave. 


m 


x 


E 


horizontal  momentum. 


|V$!Z  dz 


dx 


energy. 


(1.10) 


To  derive  the  KdV  equation  from  (1.9),  we  must  make  more 
assumptions,  to  restrict  the  possible  solutions.  For  this  restricted 
class  of  solutions,  we  may  replace  (1.9)  by  a  simpler  problem,  whose 
solutions  will  approximate  some  of  the  solutions  of  (1.9).  In 
particular,  the  KdV  equation  follows  by  making  the  following 
assumptions. 

i)  Small  amplitude  waves.  If  n  represents  a  "typical"  wave 
amplitude. 


e  *  n/h  <  1  .  (1.11a) 

11)  Long  waves.  If  k  is  a  typical  horizontal  wavenumber, 

(kh)2  «  1  .  (1.11b) 


iii)  These  two  effects  approximately  balance, 

(kh)2  ■  0(c)  .  (1.11c) 

Two  remarks  are  in  order  here.  The  first  Is  that  because  (1.9) 
uniquely  determines  its  solution  its  not  obvious  that  any  additional 
assumptions  are  permitted.  We  must  verify  a  posteriori  that  our  final 
solution  Is  consistent  with  the  assumptions  that  led  to  It.  Second, 
(kh)  and  c  are  not  well-defined.  We  may  use  them  as  convenient 
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computational  tools  only  because  they  eventually  drop  out  of  the 
analysis. 

Consistent  with  (1.11),  we  may  scale  the  independent  variables  as 
z*  *  z/h  ,  x*  =  /T  x/h ,  t*  *  /F  ct/h ,  (1.12a) 

where  c2  =  gh  (the  only  speed  available  in  the  linearized,  long-wave 
limit).  We  also  introduce  a  slow  time, 

t  =  et*  ,  (1.12b) 


so  that 

Tt  • 

Rayleigh  (1876)  noted  that  if  <p  is  analytic  at  z  =  -  h  ,  then  we 
may  expand  <p  in  a  power  series  there.  The  final  result  is  that  the 
solution  of  (1 .9a,b)  is 


[-£(1 


(l  *  z*)2]n  /_*_¥" 
(2n) !  \  3x*/ 


$0(x*,  t*.  t) 


(M3) 


If  <t>  is  analytic  at  z  ■  -  h,  this  series  is  convergent.  If  all  of 
the  derivatives  of  $0  are  bounded,  it  is  also  asymptotic  (in  e). 
Next  we  expand  the  unknowns: 

n(x,t;c)  *  h  [e  m(x*,  t*,  t)  +  c2  02  ♦  •••]  , 


(1.14) 


u 


z 


h 


3x  z  * -h 


e/gK  [u,(x\t*,T)  +  e2  u2  +  ...]  . 


It  follows  that  at  the  free  surface,  z«n 
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Now  we  simply  substitute  this  into  (1.9c)  and  the  tangential  derivative 
of  ( 1 . 9d ) ,  and  collect  terms.  The  result  at  leading  order  is: 


3m 

at* 


+ 


3Ui 

3x* 


0  , 


3ui  .  3m 

at*  lx* 


0  , 


(1.16) 


The  unique  solution  is: 

m (x*.  t*,  x )  *  f(r,  t)  +  g(t,  t)  , 

Uj(x*,  t\x)  -  f(r,  t)  -  gU.t)  , 


(1.17) 


where  r  *  x*-t*,  t  *  x*  +  t*.  The  initial  data  for  (n,  u)  define 
(f,  g),  which  inherit  their  boundedness,  smoothness,  etc.  At  this 
order,  the  solution  consists  of  a  left-  and  a  right- traveling  wave. 
There  is  no  interaction  and  no  evolution,  so  all  waves  are  permanent 
waves. 

At  the  next  order,  we  obtain 


+ 

at* 

3U;  + 

3x* 

Iru.  + 

3t 

3 

3x* 

r  (ui  hi) 

1 

6 

3 3  U  i 

3x*5 

»SLl  ♦ 

iDl  + 

lax  + 

Ul 

3ui  _  J 

3 

3u, 

at* 

3x* 

3t 

3x*  '  2 

!  3x* 

'2  3t* 

where  (ui,m)  are  defined  by  (1.17).  The  next  step  is  more 
transparent  If  we  eliminate  u2(or  n2)  and  write  the  result  In 
characteristic  coordinates: 


Similarly, 
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u2 


1 

4 


(f2-g2) 


+  \  (frr-gu>  +  -  92  • 


(1.22b) 


The  conditions  (1.21)  guarantee  that  the  left-  and  right-going 
waves  do  not  affect  each  other  long  enough  to  interact  strongly  on 
this  time-scale  [t  =  0(1)].  However,  each  of  the  two  wave  trains  does 
interact  with  itself  for  a  long  time,  and  the  two  KdV  equations 
govern  the  evolution  of  each  wave  train  on  this  longer  time-scale. 

This  is  the  main  point  of  the  derivation.  The  KdV  equation  governs 
the  evolution  on  a  slow  time-scale  of  a  small  amplitude  wave  that 
satisfied  the  linear  wave  equation  on  a  fast  time-scale. 

To  stop  at  0(e2 )  ,  one  may  set  f2  =  g2  *  0,  and  (1.22)  gives 
the  second  order  corrections  of  the  unknowns.  To  go  on  to  third  order, 
it  is  necessary  to  use  f2  and  g2  to  eliminate  secular  terms  at  the 
next  order.  In  principle,  the  expansion  can  be  carried  to  arbitrarily 
high  order,  although  this  is  rarely  done  in  practice. 

We  will  discuss  how  to  solve  the  KdV  equation  in  the  second 
lecture.  For  the  moment,  we  consider  only  some  of  the  implications  of 
this  derivation.  The  first  is  that  both  the  wave  equation  (1.16)  and 
the  KdV  equation  (1.20)  are  e-independent.  The  expansion  is  intended 
to  be  valid  in  the  limit  e  -*•  0  ,  and  we  may  take  this  limit  without 
emasculating  the  governing  equations.  That  is  indicative  that  the  KdV 
equation  is  a  true  asymptotic  equation. 

A  second  consequence  of  this  derivation  concerns  the  conservation 
laws.  It  follows  from  (1.20)  that  if  f  vanishes  rapidly  as  |r|  -*•  «  , 
or  is  periodic,  then 

I,  •  ( f  dr  ,  I2  ■  f  f2  dr  ,  I,  *  f[f’  ‘  3(fr)i:]dr 


v* 


% 
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all  are  conserved.  In  fact,  the  KdV  equation  has  an  infinite  number 
of  conservation  laws.  These  may  be  reconciled  with  the  three  physical 
conservation  laws  in  the  following  way.  We  obtained  (1.2C)  by 
expanding  the  deoendent  variables,  as  in  (1.14).  If  we  also  expand 
the  mass  integral  (for  example)  in  powers  of  c,  each  coefficient  in 
that  expansion  must  also  be  conserved.  The  infinity  of  conserved 
quantities  for  KdV  are  related  to  these  coefficients.  For  example,  we 
may  easily  verify  that  for  the  right-going  waves  (i.e.,  g  =  0),  to 
0  ( £ 2 ) , 

M  =  ^  [El!  ♦  0(£3)]  , 

mx  =  ah2  ,'gh/c  [e  1 1  +  c2  I2  +  0(c3)]  , 

(1.24) 

KE  =  [c2  1*  ♦  c3I3  ♦  0(s4)]  , 

2/e 

PE  -  WU  *  0(6*)]  , 

E  =  KE  +  PE  . 


1 
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Part  II 

Introduction  to  the  Inverse  Scattering  Transform 


The  first  lecture  was  devoted  to  the  physical  derivation  of 
equations  that  admit  solitons.  We  saw  that  even  though  the  equations 
themselves  are  fully  nonlinear,  they  typically  arise  in  physical 
problems  by  eliminating  secular  terms  in  a  weakly  nonlinear  theory. 

This  lecture  is  devoted  to  the  method  of  solution  of  these 
special  equations.  This  method  goes  under  a  variety  of  names, 
including  the  inverse  scattering  transform.  It  turns  out  that 
this  method  can  be  viewed  as  a  general ization  of  Fourier  analysis  to 
certain  nonlinear  problems.  It  provides  the  exact  solution  to 
certain  nonlinear  evolution  equations,  just  as  the  Fourier  transform 
does  for  certain  linear  evolution  equations. 

The  outline  of  this  lecture  is  first  to  review  briefly  the  method 
of  Fourier  transforms  for  linear  problems,  then  to  sketch  how  1ST 
works  for  certain  nonlinear  problems,  and  to  show  that  it  is  a 
generalization  of  Fourier  analysis.  In  the  process  we  can  say 
something  about  the  solutions  of  these  equations,  and  the  class  of 
equations  to  which  the  method  applies.  This  approach  is  essentially 
that  of  Ablowitz,  Kaup,  Newell  and  Segur  (1974). 


1. 

Linear  Evolution  Equations 

Let  us  consider  three  examples  of  linear  problems  on  the  infinite 
interval  (-«  <  x  <  ») 

(a)  ut  +  uxxx  =  0  ,  (2.1a) 

(b)  i  ut  +  uxx  =  0  ,  (2.1b) 

(c)  u-j-j-  -  u^  +  u  ■  0  ,  or 
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uxt  =  u  .  (2.1c) 

Part  (a)  may  be  called  the  Airy  equation;  it  arises  in  certain  problems 
in  optics.  Part  (b)  is  the  time-dependent  Schrodinger  equation,  with 
no  potential.  Part  (c)  is  the  Klein-Gordon  equation,  written  both  in 
laboratory  and  in  characteristic  coordinates.  In  all  three  cases,  the 
equation  holds  for  all  real  x,  initial  data  also  must  be  given  at 
t  =  0,  and  we  will  require  that  the  initial  data  be  smooth  and  decay 
rapidly  as  |x|  •*  »  . 

These  equations  all  can  be  solved  by  Fourier  transform  methods. 

The  first  step  in  that  approach  is  to  map  the  initial  data  into  its 
Fourier  transform: 


GOO 


/ 

-oo 


u(x,  0)  exp(-ikx)  dx  , 


(2.2) 


As  t  changes,  u(x,  t)  evolves  according  to  a  partial  differential 
equation,  but  u(k,  t)  satisfies  an  ordinary  differential  equation. 

(This  is  precisely  the  advantage  of  Fourier  transforms.)  The  equation 
is  so  simple  that  we  often  skip  that  step  and  simply  look  for 
solutions  in  the  form 


u(x,  t) 


2tt 


/> 


(k)  exp(ikx  -  iwt)  dk 


(2.3) 


Substituting  (2.3)  into  the  (linear)  evolution  equation  yields  the 
(linear)  dispersion  relation,  u(k)  .  In  particular,  for  our  example 
problems. 


(a)  w  s  -  k3  , 

(b)  w  *  k2  ,  (2.4) 

Ui 


(c) 


«  1/k  . 
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If  the  problem  has  a  dispersion  relation,  it  contains  all  of  the 
information  that  was  in  the  original  partial  differential  equation. 
For  example, 

(i)  for  a  first  order  (in  time)  equation,  if  m(k)  is  real  for  real 
k,  then  the  original  problem  has  an  "energy"  integral  that  is 
conserved.  In  our  example  problems, 


oo 


(ii)  If  (d2w/dk2)  f  0,  then  the  problem  is  "dispersive".  Each  wave¬ 
number  travels  with  its  own  speed,  and  the  waves  sort  themselves  out 
in  time.  [These  concepts  are  discussed  in  detail  by  Stoker  (1957) 
and  by  Whitham  (1974),  among  others.) 

The  net  effect  is  that  in  the  long  time  limit,  the  solutions  of  each 
of  these  equations  have  rather  characteristic  features  determined 
largely  by  the  group  velocity,  (du/dk)  .  For  example,  a  typical 
solution  of  (2.1a)  for  large  times  is  shown  in  Figure  2.  The  waves 
are  spreading  slowly  to  the  left,  because  (dw/dk)  <  0.  Energy  is 
conserved,  because  oj(k)  is  real,  so  as  the  waves  spread  out,  the 
amplitude  tends  to  zero  [for  (2.1a),  as  t'123  or  faster]. 

Thus,  the  solutions  of  these  conservative,  dispersive  problems 
are  comparatively  simple  for  large  times.  They  are  characterized  by 
the  waves  dispersing  over  larger  and  larger  regions  of  space,  with 
the  amplitudes  decaying  as  required  by  energy  conservation.  The 
method  to  obtain  these  solutions  depends  on  the  existence  of  two 
functions:  u(k)  represents  the  initial  data;  and  w(k)  represents 
the  evolution  equation.  This  method  of  solution  (Fourier  transforms) 
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may  be  represented  schematically  as  follows. 

u(x,0)  - -  u(k.O)  u(k,  t)  - *  u(x.t)  .  (2.6) 


2. 

Nonlinear  Evolution  Equations 

Consider  next  three  nonlinear  generalizations  of  the  equations 
in  (2.1): 

(a)  the  Korteweg-deVries  equation, 

ut  +  6uux  +  uxxx  =  0  ••  (2‘7> 

(b)  the  nonlinear  Schrodinger  equation, 

iut  +  uxx  +  2 1 u | 2  u  =  0  ;  (2.8) 

(c)  the  sine-Gardon  equation, 

uxt  *  sinu  .  (2.9) 

Why  these  examples?  First,  each  of  these  equations  linearizes  to  one 
of  the  linear  equations  in  (2.1).  Second,  they  have  the  same  energy 
integrals  as  their  linear  counterparts.  Third,  each  arises  as  a 
model  of  some  aspect  of  ocean  waves.  Fourth,  they  all  possess 
solitary  wave  solutions.  For  KdV, 

u(x,  t)  *  2k2  sech2  [<  (x  -  x0  -  4k2  t)]  ,  (2.10) 

where  k  and  xQ  are  free  constants;  see  Figure  3.  Finally,  and  most  Fg. 

importantly,  each  of  these  equations  can  be  solved  by  the  inverse 
scattering  transform  (1ST). 

Understanding  how  1ST  works  requires  some  knowledge  of  direct 
and  inverse  scattering  theory.  Here  is  an  extremely  superficial 
introduction  for  the  SchrSdinger  scattering  problem.  For  a  very 
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thorough  treatment,  see  Deift  and  Trubowitz  (1979).  Let  q(x)  be 
given,  satisfying 


J  (1  +  x2 )  |q|  dx  <  «  . 


(2.11) 


The  direct  scattering  problem  is  to  find  pairs  [X,iii(x;X)]  such  that 

ikxx  +  +  <l(x)]  =  0  ,  -oo  <  x  <  «  .  (2.12) 

Here  X  is  a  real  number,  an  "eigenvalue",  and  ip  is  required  to  be  a 
bounded  function  satisfying  certain  boundary  conditions.  For  X  >  0, 
we  may  set  X  *  k2  ,  and  require  that 

ip  ~  a (k )  exp(-lkx)  as  x  ♦  -°°  , 

(2.13) 

\p  ~  exp(-ikx)  ♦  b(k)  exp(ikx)  as  x  -  +°°  . 


These  are  related  by  | a | 2  +  | b | 2  =  1  (a  Wronskian  relation).  [In 
the  scattering  context,  these  solutions  are  called  "radiation",  a(k) 
is  called  the  "transmission  coefficient",  and  b(k)  the  "reflection 
coefficient".]  For  X  <  0,  there  are  only  a  finite  number  of 
discrete  eigenvalues  ["bound  states"].  We  may  set  X  *  -  <2  ,  and 
normalize  J  I I z  dx  *  1  •  The  boundary  conditions  are: 

-  dn  exp(<nx)  as  x  -►  -»  , 


Y  ~  cn  exp(- <n x)  as  x  -*•  +»  . 


(2.14) 


A  given  "potential",  q(x)  ,  generates  certain  scattering  data 
[a(k),  b(k);  cn,  <n].  These  may  be  collected  into  a  single  function, 

f  N 

B(x)  “  j;  J  b(k)  exp(ikx)  dk  +  exp(-<nx)  . 

n  *  1 


(2.15) 
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The  direct  scattering  problem  is  to  find  B(x)  for  a  given  q(x)  . 

The  inverse  scattering  problem  is  given  B(x)  ,  to  find  the 
potential,  q(x)  ,  that  generated  it.  This  is  an  interesting 
mathematical  question  that  was  solved  in  a  slightly  different  form 
in  a  famous  paper  by  Gel'fand  and  Levitan  (1951).  Skipping  all 
intermediate  details,  the  final  result  is  that  one  must  solve  a 
linear  integral  equation: 

OO 

K(x,  y)  ♦  B(x  +  y)  ♦  J  K(x,  z)  B(2  +  y)  dz  =  0  ,  y  >  x  . 

x 

(2.16) 

(An  equation  of  this  form  now  is  called  a  Gel ‘fand-Levi tan -Harchenko 
equation.)  Then  the  solution  of  the  inverse  scattering  problem  is 

q(x)  *  2  ■£-  K(x,  x)  .  (2.17) 

Thus,  we  may  represent  the  direct  and  inverse  scattering  problems  by 
q(x)  -  B(x  +  y )  -  K(x,y)  -  !<(x,  x)  -  q(x)  .  (2.18) 

After  that  rather  long  detour,  let  us  return  to  nonlinear 
evolution  equations,  and  come  to  the  main  point.  In  1967,  Gardner, 
Greene,  Kruskal  and  Miura  made  a  remarkable  discovery  about  the  KdV 
equation,  (2.7).  Denote  its  Initial  data  by  u(x, 0)  and  consider 

4'xx  +  [A*u(x,0)3*  *  0.  (2.19) 

In  this  way,  u(x, 0)  is  mapped  into  scattering  data,  summarized  by 
B(2x,  0)  .  As  t  changes,  u(x,  t)  evolves  according  to  (2.7),  and  of 
course,  the  scattering  data  change  as  well.  The  remarkable  fact  is 
that  if  u  satisfies  (2.7),  then 


b^a  *  o  . 


(2.20) 
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In  other  words,  u(x,  t)  evolves  through  a  family  of  potentials,  all 
of  which  have  exactly  the  same  spectrum!  [For  this  reason, 

Calogero  and  Oegasperis  (1976)  prefer  to  let  1ST  abbreviate  the 
"Iso-Spectral  Transform".]  The  rest  of  the  scattering  data  also 
evolves  simply: 


.  o. 

xr  b(k.  t)  -  8  Ik*  b  , 


_3_ 

at 


c„(t) 


4<n  cn 


(2.21) 


so  that 

B(2x,  t)  *  jjj  f  b(k)  exp [2ikx  +  81kJt]  dk 


+ 


E 


n  *  1 


c£  exp  [8*1  t  -  2<„  x]  . 


(2.22) 


Therefore  u(x,  t)  satisfies  a  nonlinear  evolution  equation  (KdV), 
but  B(2x,  t)  satisfies  a  linear  evolution  equation  whose  dispersion 
relation  Is  that  of  the  linearized  KdV  equation, 

u>(2k)  •  -8k1  .  (2.23) 

But  now  we  know  B(2x,  t)  for  any  t ,  and  we  may  reconstruct  the  new 
potential,  u(x,  t)  ,  via,  inverse  scattering.  To  sumnarlze,  we  have 
a  method  to  solve  the  KdV  equation: 

u(x,  0)  —  B(x  +y;  0)  B(x*y;  t)  —  K(x,yj  t)  -*  u(x,  t)  . 


(2.2*) 
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In  fact,  this  scheme  represents  the  method  of  solution  for  all 
problems  solvable  by  1ST.  The  steps  are: 

i)  Map  the  initial  data  of  the  nonlinear  evolution  equation  into 
the  scattering  data,  using  the  linear  scattering  problem; 

ii)  Let  the  scattering  data  evolve,  in  accord  with  the  dispersion 
relation  of  the  linearized  evolution  equation; 

iii)  Reconstruct  the  solution  of  the  nonlinear  evolution  equation 
at  a  later  time  by  solving  a  linear  integral  equation. 

Notice  that  each  step  in  this  method  is  linear,  and  that  the  whole 
procedure  parallels  the  method  of  Fourier  transforms  for  linear 
problems.  In  this  sense,  1ST  is  a  generalization  of  Fourier  analysis 
tc  certain  nonlinear  problems. 

What  about  the  solutions  of  an  equation  solvable  by  1ST  on 
-»  <  x  <  «-  ?  The  scattering  data  consists  of  a  discrete  spectrum 
(bound  states)  and  a  continuous  spectrum  (radiation),  and  these 
represent  different  kinds  of  solutions.  Each  discrete  eigenvalue 
represents  one  solitary  wave  (or  “soli ton",  since  they  are  no 
longer  required  to  be  solitary): 

B  ^  Cp  exp(-Knx)  — =*  u  -v  2**  sech2(<n(x  -  x„  -  4k*  t)}  , 


(2.25) 

where  xn  *  xn(cn)  •  The  permanence  of  these  waves  is  insured  by  the 
fact  that  the  eigenvalues  (-<*)  are  time- independent.  Notice  that 
each  of  these  waves  has  positive  speed,  and  that  the  bigger  waves 
move  faster. 

The  continuous  spectrum  requires  a  little  more  care,  but  in  a 
crude  sense  it  represents  a  part  of  the  solution  that  behaves  almost 
as  if  the  problem  were  linear.  More  precisely,  it  represents  a 
dispersive  wave-train  in  which  the  linearized  group-velocity  (which 
is  negative  for  KdV)  plays  an  Important  role.  The  net  result  Is  that 
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the  long-time  solution  of  KdV  is  comparatively  simple,  as  shown 

in  Figure  4.  It  may  be  worth  emphasizing  that  even  though  the  Fg.4 

equation  is  fully  nonlinear,  its  long  time  behavior  may  be  predicted 
to  any  desired  accuracy.  The  problem  is  stable,  and  there  are  no 
chaotic  solutions. 


3. 

Generalizations 

The  reader  may  have  the  impression  by  now  that  1ST  is  a  miracle 
that  works.  To  a  certain  extent  that  is  true,  but  that  is  no  reason 
not  to  use  it  wherever  possible.  The  question  is,  "To  what  problems 
does  1ST  apply?"  Next  we  show  that  there  are  infinitely  many 
problems  solvable  by  1ST.  (So  even  though  it  is  a  miracle,  it  may  not 
be  uncommon.)  Let  us  consider  a  different  scattering  problem  on 

-oo  <  X  <  oo  ; 


3X  Vi  +  i  c  Vi  -  q  V2  , 
3XV2  -  i  qVj  *  rVi  , 


(2.26) 


with  q,  r  -►  0  as  |x|  -►  ®.  This  problem  was  first  analyzed  by 

Zakharov  and  Shabat  (1972)  for  r  =  - q*  (*  denotes  complex  conjugate). 

Note  that  if  r  *  1  ,  (2.26)  can  be  reduced  to  (2.12).  Note  further 

that  if  r  *  0  and  Vj  exp(icx)  ■»  1  as  x  *  -» ,  then 

lim  [V,  exp(i  ;x)]  Is  just  the  Fourier  transform  of  q(x)  .  In 
x -*■+»> 

this  sense  (2.26)  is  already  a  generalization  of  Fourier  transforms. 

To  construct  an  1ST  out  of  (2.26),  we  need  [q(x,  t),  r(x,  t)]  to 
evolve  so  that  3t c  ■  0 ;  l.e.,  we  will  force  the  eigenvalue  to  be 
time-independent.  To  do  this,  we  allow  the  eigenfunctions  to  evolve 
according  to  linear  equations: 
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where 

that 


* 

AVi  +  BV2  , 

(2.27) 

3tV2  = 

c  V,  -  av2  , 

A  *  A(q,r;  c)  ,  etc.  Compatibility  of  (2.26)  and 

(2.27)  requires 

(V)xt  *  (V)tx  •  Demanding 

3t  c  =  0  yields: 

Ax 

=  q  C  -  r  B  , 

Bx  +  2i  z  B 

■  qt  -  2Aq  , 

(2.28) 

Ct  -  21  c  C 

*  rt  +  2A  r  , 

If  (A,  B,  C)  satisfy  these  coupled  ordinary  differential  equations 
(in  x),  then  the  eigenvalue  is  constant,  and  we  can  construct  an  1ST 
based  on  (2.26).  Boundary  conditions  for  (2.28)  are  obtained  by 
comparing  (2.26),  (2.27)  as  |x|  -*■  « .  Because  (q,  r)  -0,  we  know 
that  Vj  n,  Ci  exp(-isx)  ,  V2  ■x-  C2  exp(icx)  from  (2.26).  To  assure 
compatibility  with  (2.27),  we  require 

A  -  A0U)  ,  B,C  -  0  as  |x|  -  »  .  (2.29) 

This  gives  six  boundary  conditions  for  (2.28),  which  is  now  over- 
determined.  Therefore  (2.28),  (2.29)  will  have  no  solution  except  in 
special  cases. 

Given  (2.26),  the  choice  of  nonlinear  evolution  equations  is 
determined  by  A  (;)  .  Here  we  use  the  linearized  dispersion  relation: 

»(2C)  ■  2iA0(-5)  .  (2.30) 

For  example,  if  A0(c)  -  -  21  c2  .  then  (2.28)  has  a  solution  only  if 
1qt  +  qxx  -  2q2r  -  0  , 

i  rt  -  rxx  +  2  q  rJ  -  0  . 


(2.31) 
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lf  r  =  -  q*  ,  each  of  these  reduce  to  the  nonlinear  Schrodinger 
equation,  (2.8). 

Theorem  (Ablowitz,  Kaup,  Newell,  and  Segur,  1974) 

Let  u(k)  be  a  ratio  of  entire  functions,  and  be  real 
for  real  k  .  Then  (2.26)  generates  a  nonlinear  evolution 
equation  that  is  solvable  by  1ST,  and  whose  linearized 
dispersion  relation  is  w(k)  .  If  w(k)  also  is  an  odd 
function  of  k,  then  (2.12)  generates  a  different  nonlinear 
equation  solvable  by  1ST.  Its  linearized  dispersion  relation 
also  is  u(k)  . 

If  your  objective  is  to  generate  nonlinear  equations  solvable  by 
1ST,  you  need  two  ingredients: 

i)  A  scattering  problem;  and, 

ii)  A  linearized  dispersion  relation. 

Each  such  pair  generates  one  nonlinear  problem  solvable  by  1ST.  If 
your  objective  is  to  solve  a  particular  equation,  this  is  almost  no 
help  at  all.  In  my  fifth  lecture,  I  will  discuss  how  to  determine 
whether  a  given  equation  can  be  solved  by  1ST. 


r 


i 
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Part  III 

More  Inverse  Scattering  on  the  Infinite  Interval 


In  the  second  lecture  we  saw  that  the  inverse  scattering 
transform  works,  and  that  it  may  be  regarded  as  a  generalization  of 
the  Fourier  transform  to  certain  nonlinear  problems.  For  the  KdV 
equation  (2.7),  1ST  may  be  represented  schematically  as  follows: 


direct 

scattering 

u(x,  0)  - ►  {a ( k,  0),  b ( k ,  0),  ,  cn(0) } 


M  =  o 

1  ab  =  , 

at 

b  at  ‘  ' 

1  ^ 

ii 

D 

^  3cn  _ 

at 

cn  3t 

(3.1) 


u(x,  t)  i - (a(k,  t),  b(k,  t),  k  ,  cn( t) } 

inverse 

scattering 


In  this  lecture  we  will  examine  1ST  on  the  infinite  line  in  more 
detail.  The  extra  information  may  help  to  explain  why  1ST  works  and 
what  sort  of  solutions  it  admits. 


1. 

Hamiltonian  Mechanics 

Zakharov  and  Faddeev  (1971)  pioneered  a  description  of  1ST  as  a 
canonical  transformation  of  a  Hamiltonian  system  to  action-angle 
variables.  This  description  is  an  alternative  to  that  of  1ST  as 
Fourier  analysis  for  nonlinear  problems.  Both  are  legitimate;  which 
is  preferable  is  a  matter  of  taste. 

The  reader  may  recall  that  Hamiltonian  mechanics  is  simply  a 
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variational  description  of  certain  dynamical  systems.  (Basic 
references  here  are  Arnold,  1978;  or  Goldstein,  1950.)  In  this 
formulation,  one  identifies  generalized  "coordinates"  (q)  and 
"momenta"  (p),  which  describe  completely  the  state  of  the  system.  If 
the  problem  has  N  degrees  of  freedom,  then  p  and  q  are  each 
N-dimensional  vectors.  For  the  partial  differential  equations  of 
interest  here,  p  and  q  are  infinite-dimensional.  One  introduces  a 
Hamiltonian,  H(p,q),  which  must  have  the  property  that  Hamilton's 
equations, 


q 


(3.2) 


are  equivalent  to  the  equations  of  motion  of  the  system.  Here 
(")  =  3t(  )  >  and  the  derivatives  of  H  in  (3.2)  are  functional 
derivatives.  In  general  H  may  depend  on  time,  but  for  any  of  the 
conservative  systems  under  consideration  it  does  not.  Any  system 
that  has  such  a  variational  formulation  is  said  to  be  Hamiltonian. 

It  happens  that  problems  solvable  by  1ST  are  Hamiltonian.  As 
examples,  consider 


H  '  -  j *Pqx  +  P2qx  -  px  qxx  1  dx- 

h2  *  - 1  f  I  qx  px  +  p‘  q2  ( dx . 


(3.3) 


Hj  yields  two  evolution  equations,  which  admit  the  identification 
p  =  qx  .  Each  equation  reduces  to  KdV  under  this  identification. 
Similarly,  H2  gives  the  nonlinear  Schrodinger  equation  if  we 
identify  p  =  t  q*  (complex  conjugate). 

Canonical  transformations  play  an  important  role  in  Hamiltonian 
mechanics.  Roughly  speaking,  a  canonical  transformation  is  simply  a 
change  of  variables. 


(P.  q) 


(P,Q) 


(3.4) 


that  preserves  the  volume  of  the  phase  space;  (this  corresponds  to  a 
mapping  that  is  1-1  and  onto).  One  may  check  this  property  either  by 
means  of  Poisson  brackets  (old-f3Shior.ed)  or  a  symplectic  form 
(new-fangled).  Necessarily,  a  canonical  transformation  does  not 
affect  the  form  of  Hamilton's  equations: 

H(P,Q)  :  0  -  §  ,  P  =  -  ■  (3.5) 

Out  of  all  possible  canonical  transformations,  an  especially 
desirable  one  is  one  in  which  the  new  Hamiltonian  is  independent  of  Q. 
These  are  called  action  (P)  and  angle  (Q)  variables.  Obviously,  if 
H  =  H(P),  then  from  (3.5), 

P  -  0  ,  Q  =  =  constant  ,  (3.6) 

5p 

so  that  integration  of  the  equations  is  as  simple  as  possible. 

Another  way  to  say  this  is  that  if  a  Hamiltonian  system  has  action- 
angle  variables,  then  its  motion  is  basically  very  simple,  when 
viewed  appropriately.  Unfortunately,  there  is  no  general  method  known 
to  determine  whether  a  system  has  action-angle  variables,  or  to  find 
them  if  they  exist. 

The  valuable  insight  of  Zakharov  and  Faddeev  (1971)  was  that  (3.1), 
the  equations  that  describe  the  evolution  of  the  scattering  data  for 
KdV,was  in  the  form  of  (3.6). 

Theorem  (Zakharov  and  Faddeev,  1971) 

Let  u(x,  t)  represent  a  KdV  solution  on  -»  <  x  <  ». 

The  1ST  mapping 

P(k)  «  k in | a ( k ) ] 2  ,  Q(k)  =  Im  in  b(k) 

u(x,  t)  - - 

=  <n  ’  =  2,0  cn 

(3.7) 

is  a  canonical  transformation  to  action-angle  variables. 


I 
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From  this  standpoint,  it  is  not  surprising  that  the  KdV  equation 
has  an  infinite  set  of  constants  of  the  motion.  They  are  a 
representation  of  the  (time-independent)  action  variables.  Moreover, 
given  that  the  equation  on  -=>  <  x  <  =>  has  an  asymptotic  state,  it  is 
not  surprising  that  the  asymptotic  behavior  is  relatively  simple.  The 
existence  of  action-angle  variables  means  that  motion  is  basically 
simple  (when  viewed  appropriately);  the  simple  asymptotic  behavior 
reveals  the  basically  simple  motion. 


2. 

Scattering  Theory 

Another  apparent  miracle  in  1ST  is  the  fact  that  scattering 
theory  works  so  well.  Even  after  accepting  that  the  potential  in  the 
scattering  problem  [e.g.,  u(x)  in  (2.12)]  is  determined  by  the 
scattering  data,  one  is  still  surprised  at  the  simplicity  of  the 
inversion  procedure.  In  fact,  scattering  theory  works  as  well  as  it 
does  because  it  uses  the  powerful  theory  of  analytic  functions  of  a 
complex  variable. 

We  may  illustrate  this  close  connection  to  the  theory  of 
analytic  functions  by  considering  a  scattering  problem  due  to 
Zakharov  and  Shabat  (1972),  and  generalized  by  Ablowitz,  Kaup,  Newell 
and  Segur  (1974): 


3xVj  +  i  sVi  =  q(x)  V2  , 
3XV2  -  i  ?V2  =  r(x)  V,  . 


(3.8) 


Here  -»  <  x  <  ®  ,  and  we  assume  that  q(x),r(x)  vanish  rapidly  as 
| x |  •*  ®  .  This  scattering  problem  is  appropriate  for  the  nonlinear 
Schrodinger  equation,  the  sine-Gordon  equation  and  infinitely  many 
other  problems.  However,  as  time-dependence  is  not  germane  to 
scattering  theory  £er  se ,  we  will  hold  time  fixed  for  the  current 
discussion.  Thus  (q,  r)  may  be  considered  known  functions  that  are 
absolutely  integrable. 
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(3.9) 


For  real  C,  $  and  41  are  linearly  independent  for  all  real  x,  as 
are  4  and  4.  Because  (3.8)  has  only  two  linearly  independent 
solutions,  we  may  define  (a(c)  ,  a(c)  ,  b(c)  ,  b(c)}  by 


4  =  a4  +  b4  , 
4  =  b?  +  a4  . 


(3.10) 


These  are  related  by  a  Wronskian  relation,  aa  +  bb  =  1  ,  From  (3.9) 
and  (3.10), 


The  set  {a,  a,  b,  b}  for  real  c  makes  up  part  of  the  scattering  data. 
These  functions  may  be  viewed  as  representing  the  asymptotic  behavior 
of  certain  solutions  of  (3.8)  as  in  (3.11),  or  simply  as  Wronskians 
of  solutions  from  (3.10). 

Once  these  functions  have  been  defined  for  real  c.  they  may  be 
extended  into  the  complex  i-plane.  Ablowitz,  Kaup,  Newell  and 


Segur  (1974)  proved  that: 


i)  [+>( x,  ;)  exp(i;x)]  and  [>(x,  t)  exp(-i;x)]  are  analytic 
functions  of  ;  for  all  real  x  if  Im(;)  >  0  ; 

ii)  [4  exp(-i;x)]  and  [?  exp{i;x)]  are  analytic  in  5  if 
Im(;)  <  0  ; 

iii)  a(;)  -  lim  41  exp(Ux)  is  analytic  for  Im(?)  >  0,  and 
x  -=> 

a(;)  is  analytic  for  Im(4)  <  0; 

iv)  the  discrete  eigenvalues  of  (3.8)  are  zeros  of  ak)  or 
ak)  in  their  regions  of  analyticity; 

v)  as  |;|  +  » ,  Imk)  >  0  ,  ak)  -  1  +  Ok'1 )  , 

02  exp(kx)  —  -  [r(x)/2i?]  +  0(c~2 )  , 

i|n  exp(-kx)  — ►  [q(x)/2ic]  +  Ok'2)  ,  with  similar  results 
for  Imk)  <  0  . 

Most  of  the  results  of  1ST  are  a  consequence  of  these  relations 
and  the  time-dependence  of  the  scattering  data.  In  particular,  one 
finds  that  ak)  and  ak)  are  time-independent.  Then  it  follows 
from  (iv)  that  the  discrete  eigenvalues  are  time-independent  as  well 
(the  iso-spectral  property).  It  follows  from  (v)  that  in  a( c )  has 
an  asymptotic  expansion  as  |;|  ♦  “  for  Im(c)  >  0: 

00 

inak)  =  2  In/(2ic)n  •  (3.12) 

n  =  1 


Because  a (?)  is  time-independent,  the  coefficients  in  this  expansion 
must  be  time-independent  as  well.  These  are  the  infinite  set  of 
conserved  integrals  of  any  of  the  evolution  equations  solved  by  (3.8). 
It  also  follows  from  (v)  that  if  we  can  reconstruct  $2  exp(iix)  and 
4>i  exp(-icx)  from  the  scattering  data,  then  q(x)  and  r(x)  may  be 
obtained  by  taking  a  limit.  We  show  next  that  inverse  scattering 
theory  follows  just  that  strategy. 
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To  obtain  the  linear  integral  equations  that  are  the  heart  of 
inverse  scattering  theory,  write  (3.10a)  as 

exp(i,xj^  _  ^  exp(i;x)  +  exp(i;x)  ,  q  real  ,  (3.13) 

3  3 

For  simplicity,  let  us  assume  that  a(q)  has  no  zeros  for  Im(;)  >,  0. 
Then  [j>  exp(ijx)]/a  is  analytic  for  Im(;}  >  0  and  vanishes  as 
| c |  -*•  “  there,  while  4>exp(i;x)  has  similar  properties  for 
Im(;)  <  0  .  Now  the  problem  of  reconstructing  these  analytic  functions 
from  (3.13)  is  very  similar  to  a  famous  problem  posed  by  Hilbert 
(cf.,  Muskhelishvili,  1953,  Ch.  5).  That  problem  may  be  stated  as 
follows.  F+(z)  is  analytic  for  Im(c)  >  0  and  vanishes  as 
j z |  •*  °°  there.  F_ ( z )  is  analytic  for  Im(;)  <  0  and  vanishes 

as  j z |  -  ®  there.  They  are  related  on  the  real  axis  by: 

F+(x)  -  F_(x)  =  f(x)  on  z  =  x  ,  (3.14) 


where  f(x)  is  a  given  function.  The  "Hilbert  problem"  is  to 
construct  both  F+  and  F_  from  f(z)  .  If  f(x)  is  absolutely 
integrable,  the  solution  of  the  problem  is  given  by: 


F(z) 


_L  f  f(x)dx  . 

2iti  J  x  -  z 


F+(z)  , 
F_(z)  , 


Im(z)  >  0  , 

(3.15) 

Im(z)  <  0  . 


Comparing  (3.13)  and  (3.14),  it  is  evident  that  because  (b/a)  is 
known  for  real  c,  if  ijjexp(icx)  were  known  for  real  c,  then  (3.13) 
would  be  solved  by  a  formula  like  (3.15).  However,  because  ijiexp(itx) 
is  unknown,  we  obtain  instead  a  linear  integral  equation,  with  a 
Cauchy-type  singularity.  Another  singular  integral  equation  follows 
from  (3.10b).  The  usual  Gel 'fand-Levi tan- type  of  integral  equations 
for  (3.8)  are  essentially  the  Fourier  transform  of  these  coupled 
singular  integral  equations. 

The  original  work  by  Gel'fand  and  Levitan  (1951)  and  others  on 
inverse  scattering  theory  did  not  make  this  connection  to  the  Hilbert 
problem,  which  has  been  developed  more  recently  (e.g.,  Zakharov  and 
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Manakov,  1979).  The  advantage  of  this  approach  is  that  it  emphasizes 
the  fundamental  role  played  by  analytic  functions  in  the  theory  of 
inverse  scattering. 


3. 

Solutions  of  the  Nonlinear  Schrodinger  Equation 

The  practical  consequence  of  all  of  this  remarkable  structure  is 
that  the  solutions  of  these  special  nonlinear  evolution  equations  are 
quite  predictable,  and  can  be  computed  explicitly  (especially  for 
large  times)  once  the  initial  data  has  been  mapped  into  scattering 
data.  We  may  illustrate  this  by  focussing  our  attention  on  the 
nonlinear  Schrodinger  equation, 

1  qt  +  qxx  +  2lql2q  =  0  •  (3.16) 

which  describes  the  nonlinear  instability  of  a  packet  of  nearly 
monochromatic  water  waves  of  small  amplitude  in  one  dimension. 

The  appropriate  scattering  problem  is  (3.8),  with  r  *  -  q*  . 

The  discrete  spectrum  in  this  case  is  represented  by  N  discrete 
eigenvalues,  *  (£  +  in)j  with  nj  >  0  ,  along  with  N  parameters, 

Cj  .  The  continuous  spectrum  is  represented  by  b/a  (£)  ,  on  n  =  0  . 
[The  equation 

i qt  +  qXx  ■  2lql2q  "  0  •  (3.17) 

also  has  physical  interest.  In  this  case,  r  ■  +q*  in  (3.8),  and 
there  are  no  discrete  eigenvalues.] 

A  single  eigenvalue  with  n  >  0  corresponds  to  a  single  soliton 
solution  of  (3.16). 

q(x,  t)  *  2n  sech  {2n(x  +  4?t  +  Xq))  exp  {-  2i[£x  +  2(£2  -  n2)  t  +  <*>„]}  , 

(3.18) 


This  is  an  "envelope  soliton"  (as  opposed  to  a  KdV-type  soliton).  It 
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represents  a  one-dimensional  wave  packet  that  is  stable  with  respect 

to  one-dimensional  perturbations,  and  is  shown  in  Figure  5.  For  a  Fg.  i 

one-dimensional  packet  of  water  waves  that  are  nearly  monochromatic 

and  of  small  amplitude,  (3.16)  describes  the  Benjamin-Feir  (1967) 

instability,  and  one  could  think  of  an  envelope  soli  ton  as  a  stable 

equilibrium  state  for  that  process. 

N  discrete  eigenvalues  generate  N  solitons.  If  they  have  N 
different  speeds  (£j  /  ^  for  j  f  k),  they  separate  in  space  as 
t  •*  <=  ,  so  that  the  long-time  solution  of  (3.16)  is  simply  a  sum  of 
N  individual  solitons.  We  should  emphasize  that  solitons  are 
intrinsically  nonlinear  objects,  that  disappear  in  the  linear  limit. 

On  the  other  hand  the  "radiation",  which  corresponds  to  the 
continuous  spectrum,  is  qualitatively  quite  similar  to  the  solution 
of  the  linearized  problem.  For  comparison  we  note  that 

iyt  >  yxx  -  o  (3.19) 

has  a  family  of  self-similar  solutions,  including 

y  =  t'*5  A  exp(i  x2/4t  +  i  <(>)  .  (3.20) 

One  may  also  solve  (3.19)  as  an  initial  value  problem  on  -=>  <  x  <  » , 
and  evaluate  the  solution  in  the  long  time  limit.  The  result  is  that 
along  x/t  =  2k  (the  group  velocity),  * 

y(x,  t)  ~  t'*5  Y(k)  exp(i  x2/4t  -  1  tt/4)  ,  (3.21) 

where  9 ( k )  is  the  Fourier  transform  of  the  initial  data.  This  has 
the  form  of  a  slowly-varying  similarity  solution. 

There  is  a  nearly  Identical  similarity  solution  of  (3.16), 

q  =  t'*5  A  exp[1(x2/4t  +  2A2  in  t  +  0)]  .  (3.22) 

In  the  absence  of  solitons,  the  asymptotic  solution  of  (3.21)  that 
evolves  from  appropriate  initial  data  also  takes  the  form  of  a 
slowly-varying  similarity  solution.  Along  x/t  *  -  4£  , 

— 


X 
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q  ~  t-"5  f(x/t)  exp  {i[x2/4t  +  2f2  Unt  +  g(x/t)]}  ,  (3.23) 


where 

f2(x/t)  =  -47  *"[1  -  Ib/a(^)|2]  ,  (3.24) 

and  g(x/t)  also  is  determined  by  b/a  (£)  .  [See  Zakharov  and 
Manakov,  1976;  Ablowitz  and  Segur,  1975;  and  Segur,  1975.] 

The  general  solution  of  (3.16)  on  -®  <  x  <  ®  involves  both 
soli  tons  and  radiation.  For  KdV  these  two  components  separate  in 
space,  but  for  (3.16)  they  coexist.  In  the  long  time  limit,  the 
solution  consists  of  N  envelope  solitons  riding  on  a  sea  of 
radiation.  The  solitons  are  permanent  wave  packets,  while  the 
radiation  decays  as  t  . 
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Part  IV 

The  Korteweg-deVries  Equation 
with  Periodic  Boundary  Conditions 


To  this  point  we  have  discussed  solitons  and  1ST  only  on  the 
infinite  interval.  However,  the  KdV  equation  with  periodic  boundary 
conditions  also  has  applications  in  water  waves.  Moreover,  the 
original  discovery  of  solitons  by  Zabusky  and  Kruskal  (1965)  was 
based  on  numerical  experiments  on  the  periodic  KdV  problem. 

There  is  a  theory  of  inverse  scattering  transforms  for  the  KdV 
equation  with  periodic  boundary  conditions.  The  most  complete 
version  is  due  to  McKean  and  Trubowitz  (1976).  In  this  lecture,  we 
will  follow  the  less  general  but  simpler  version  of  Dubrovin  and 
Novikov  (1974).  In  contrast  to  the  theory  of  1ST  on  the  infinite 
interval,  however,  the  theory  for  the  periodic  problem  cannot  yet  be 
considered  a  practical  tool  for  use  in  applications.  After  presenting 
the  theory  in  its  current  form,  we  will  identify  some  practical 
questions  that  have  not  yet  been  answered  satisfactorily. 

The  first  work  on  the  periodic  KdV  problem, 


ut  +  6uux  +  uxxx  *  0  , 

u(x,  t)  =  u(x  +  L,  t)  , 


(4.1) 


was  done  by  Korteweg  and  deVries  (1895),  who  found  a  periodic, 
traveling  wave  solution  of  (4.1): 

u(x,  t)  *  2p2  k2  cn2[p(x  -  ct  +  x);  k]  +  S  .  (4.2) 

Here  cn(<pj  k)  is  the  Jacobian  elliptic  function  with  modulus  k 
(0  <  k2  S  1)  , 


c  *  6S  -  4p2(l  -  2k)  , 
PL  * 


2K(k)  , 


(4.3) 
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where  K(k)  is  the  complete  elliptic  integral  of  the  first  kind 
(cf.,  Byrd  and  Friedman,  1971).  If  we  require  ^udx  =  0,  as 
required  by  (1.21),  then 

3  ■  -2'1  [ft}-'* k’]  •  (4-" 

where  E ( k )  is  the  complete  elliptic  integral  of  the  second  kind. 

Korteweg  and  deVries  called  these  "cnoidal  waves",  by  analogy  with 

sinusoidal  waves.  A  typical  solution  is  shown  in  Figure  6,  for  Fg.  6 

k  =  4  .  These  nonlinear,  periodic,  traveling  waves  reduce  to 
infinitesimal  sinusoidal  waves  if  k  ■*  0  , 

u(x,  t)  •*  p2  k2  cos  2p(x  +  4p2  t  +  x)  ,  (4.5) 

with  pL  -  w  .  At  the  other  extreme,  if  k  •*  1  , 

u(x,  t)  -►  2p2  sech2  p(x  -  4p2  t  +  x)  .  (4.6) 

Aside  from  these  special  solutions,  however,  almost  nothing  was 
known  about  the  periodic  KdV  problem  until  the  important  numerical 
work  of  Zabusky  and  Kruskal  (1965),  who  were  motivated  by  the  earlier 
work  of  Fermi,  Pasta  and  Ulam  (1955).  In  both  of  these  studies,  the 
authors  observed  "recurrence  of  initial  states"  after  a  relatively 
short  time.  In  other  words,  rather  than  showing  any  tendency  to  an 
equi partition  of  energy  among  all  of  the  degrees  of  freedom  of  the 
system,  the  solution  of  the  initial  value  problem  almost  returned  to 
its  initial  configuration  repeatedly.  The  time  required  for  this 
recurrence  was  short  enough  to  suggest  that  only  a  few  of  the 
possible  degrees  of  freedom  of  the  system  actually  were  participating 
in  the  process  described  by  (4.1).  But  how  could  the  solution  of 
(4.1)  be  so  constrained  unless  (4.1)  Itself  carried  additional 
constraints  that  had  not  been  discovered?  These  extra  constraints 
turned  out  to  be  the  infinite  set  of  conservation  laws  of  the  KdV 
equation  (Mlura,  Gardner  and  Kruskal,  1968).  The  discovery  of  these 
conservation  laws  led  in  turn  to  the  development  of  1ST  (Gardner, 
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Greene,  Kruskal  and  Miura,  1967).  We  should  emphasize  the  historical 
importance  of  the  careful  numerical  studies  of  Zabusky  and  Kruskal 
(1965),  which  indicated  that  the  KdV  equation  possessed  additional 
mathematical  structure. 

Let  us  now  outline  the  theory  of  1ST  for  the  KdV  problem.  The 
reader  will  observe  that  some  of  this  theory  is  analogous  to  that  on 
the  infinite  interval,  but  some  of  it  is  not.  The  scattering  problem 
is: 


Vxx  +  [A  +  u(x)]V  =  0  ,  (4.7) 

just  as  it  was  on  the  infinite  interval,  but  in  this  case  u(x)  is  a 

periodic  function  of  x.  Because  scattering  theory  requires  no 
knowledge  of  the  time-dependence  of  u(x,  t)  ,  we  may  consider  (4.7) 
at  a  fixed  time. 

One  of  the  difficult  conceptual  questions  about  the  periodic 
KdV  problem  is  to  find  what  plays  the  role  of  the  scattering  data 
here,  since  there  is  no  "point  at  infinity"  where  the  solution 
simplifies.  Without  such  a  special  point,  we  simply  choose  an 
arbitrary  xQ  in  order  to  begin  the  analysis.  With  xQ  fixed  (and 
time  fixed),  we  may  identify  two  linearly  independent  solutions  of 
(4.7),  4>(x;  xQ,  X)  and  its  complex  conjugate,  $*,  by  imposing 

boundary  conditions  for  real  A  at  x  «  xQ  : 

<t>(x0;  x0,  X)  =  1  d>* ( x0 ;  Xq,  X)  =  1 

$x(x0;  xq,A)  =  ik  «  ivT  bj(x0;  xo*  =  *  ik 

(4-8) 

One  period  to  the  right  (x  •*  x  +  L),  these  two  functions  satisfy  the 
same  differential  equation  again,  so  they  must  be  a  linear  combination 
of  the  <t>  and  : 


<t>(x  ♦  L;  x0,  X) 

O 

X 

*0 

1 _ 

b(Xg,  A) 

4>(x;  x0,  A) 

$*(x  +  L;  Xg ,  X) 

b*(Xg, A) 

a*(x0.  A) 

$*(x;  x0.  A) 

(4.9) 
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The  matrix  connecting  these  two  sets  of  solutions  of  (4.7)  is  called 
the  "monodromy  matrix".  In  some  ways  it  plays  the  role  of  the 
scattering  data.  Its  coefficients  are  related  by  a  Wronskian 
relation. 


=  1 


(4.10) 


The  fact  that  the  coefficients  of  (4.6)  are  periodic  does  not 
necessarily  mean  that  its  solutions  are  periodic.  In  fact,  most  of 
them  are  not.  The  solutions  of  (4.6)  that  are  periodic  play  a 
fundamental  role  in  t  a  theory  of  inverse  scattering.  We  define  next 
the  "Bloch  eigenfunctions"  to  be  solutions  of  (4.6)  that  satisfy 


<k(xQ  »  Xq  ,  X )  -  1  , 
ili(x  +  L;  Xq ,  X)  =  u  ip ( x ;  xQ,  X)  . 


(4.11) 


Because  these  must  also  be  a  linear  combination  of  $  and  ,  one 
may  show  that  for  each  (xQ,  X)  , 

u2  -  2aru  +  1  =  0,  (4.12) 


where 


a  =  ar  +  i  a^  .  (4.13) 

Equation  (4.12)  admits  three  possibilities. 

i)  If  |ar|  >  1  ,  one  root  of  (4.12)  is  larger  than  one  in 
magnitude,  and  the  other  is  smaller.  From  (4.11),  these 
correspond  to  Bloch  eigenfunctions  that  grow  without  bound, 
either  as  x  ♦  +=»  or  as  x  •+  -»  .  These  are  said  to  be 
"unstable" . 

ii)  If  | ar |  <  1  ,  we  may  define  ar(X)  *  cosp(x)  ,  and  show  that 


u  *  exp(ip)  . 


(4.14) 


These  represent  "stable"  eigenfunctions, 
i  i  i )  I f  ar !  =  1  ,  i .  e . , 

a‘  -  1  ,  (4.15) 

then  u  *  ±1  and  the  eigenfunctions  are  either  periodic  or 
antiperiodic  functions  of  x.  Thus,  (4.15)  identifies  the 
periodic  solutions  of  (4.6). 

These  possibilities  may  be  summarized  on  a  "Floquet  diagram", 
in  which  ar  is  plotted  as  a  function  of  X  for  fixed  x0  ,  as 
shown  in  Figure  7.  The  regions  in  which  |ar|  >  1  are  called 
"unstable  bands".  Necessarily,  each  unstable  band  lies  between  two 
successive  roots  of  (4.15),  numbered  X2n  and  ^n+l  •  0ne  ma^  show 
(by  oscillation  theorems;  cf.,  Magnus  and  Winkler,  1979)  that 

X2n  >  *2n-l  ’  and  that  x2n+l  -  X2n  '  Moreover>  one  stationary 
point  of  ar  occurs  in  each  unstable  band,  and  none  occur  elsewhere. 
Unstable  bands  that  consist  of  a  single  point  (if  X2p  *  X2n+1 ^  are 
called  "degenerate". 

Following  Dubrovin  and  Novikov  (1974),  we  now  assume  that  u(x) 
generates  ar(xQ,  X)  with  only  a  finite  number  (N)  of  non-degenerate 
bands.  (Roughly,  this  corresponds  to  approximating  a  periodic 
function  with  a  finite  Fourier  series.) 

Next,  we  define  two  spectra  for  (4.6).  These  may  be  defined 
either  by  attaching  boundary  conditions  to  (4.6),  or  by  imposing 
conditions  on  the  monodromy  matrix. 

i)  The  main  spectrum  is  defined  by  (4.15)  or  by  requiring  that 
(4.6)  admit  periodic  solutions.  Points  of  this  spectrum 
(Xj,  X2,...)  occur  at  the  edges  of  the  unstable  bands.  Note 
that  periodic  solutions  of  (4.6)  are  independent  of  xQ  . 

ii)  The  auxiliary  spectrum  (-n,  >2,...)  is  defined  by 


or  by  requiring  that  there  be  a  solution  of  (4.6)  that  satisfies 


y(xQ;  xQ,  X)  =  0  «  y0(x0+L;  xQ,  X) 


It  follows  from  (4.10)  and  (4.16)  that  a£  >  1  at  each  of  these 
points,  which  must  therefore  lie  in  the  unstable  bands.  One  may  show 
that  each  unstable  band  contains  exactly  one  point  of  the  auxiliary 
spectrum.  By  assumption,  at  most  N  of  these  points  of  the  auxiliary 
spectrum  do  not  coincide  with  points  of  the  main  spectrum. 

Given  u(x)  and  xQ  ,  determination  of  (Xj)  and  (yj) 
completes  the  direct  scattering  problem.  The  inverse  mapping 
requires  knowledge  of  the  analytic  properties  of  the  monodromy  matrix, 
as  one  might  expect  by  analogy  with  the  inverse  problem  on  the 
infinite  interval.  The  final  result  is  miraculously  simple: 


2N  +  1  N 

n  +  n  . 

j  - 1  J  j  - 1  J 


Now  it  remains  only  to  find  how  the  two  spectra  depend  on  xQ  ,  and 
on  time. 

The  boundary  conditions  that  define  (Xj)  and  (yj)  suggest 
that  the  main  spectrum  should  be  independent  of  x0  ,  but  that  the 
auxiliary  spectrum  should  not  be.  This  turns  out  to  be  the  case.  One 
shows  that 


(4.18a) 


j  -  1 . N 


(4.18b) 


where  Oj  *  ±  1  and 


(4.1  Sc } 


R(-) 


2N  +  1 

n  {■•  -  -j) 
j  =  i  J 


The  x0-decendence  of  -C v  . >  is  as  frightening  as  (4.17)  was  appealing. 
One  is  faced  in  (4.13b)  with  the  integration  of  N  coupled,  nonlinear, 
ordinary  differential  equations,  even  before  time-dependence  is 
brought  into  the  picture!  Fortunately,  another  miracle  occurs,  and 
there  is  a  change  of  variables  (involving  hyperel 1 iptic  functions)  that 
permits  one  to  integrate  (4.13b)  by  quadrature.  Thus  although  (4.13b) 
looks  foreboding,  it  represents  a  straight-line  motion  when  viewed 
appropriately.  The  cost  of  acquiring  this  simple  picture  is  that 
one  must  introduce  hyperel 1 iptic  functions. 

The  time-dependence  of  the  two  spectra  parallels  their 
xQ-dependence.  For  the  KdV  equation, 


it 


(4.19a) 
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R'-(fj) 
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1,...,  N 


(4.19b) 


where  *j  =  ±  1  ,  and  R(a)  was  defined  by  (4.18c).  Again,  the 
integration  in  (4.19b)  may  be  reduced  to  quadrature  by  introducing 
hyperel 1 iptic  functions. 

What  are  the  consequences  of  introducing  hyperelliptic  functions? 
One  may  show  that  for  an  N-band  solution  of  the  KdV  equation  (i.e., 
whose  auxiliary  spectrum  contains  only  N  nondegenerate  points), 


u(x,  t) 


^  jx2  9^(hl  ,  12 


. . ,  hN)  +  constant  , 


(4.20a) 


where 


I 


1- 


- f 


ij  =  i (kj  x  -  -j  t)  ,  (4.20b) 

each  kj  and  is  constant,  and  9^  is  the  theta  function,  an 

analytic  function  that  is  periodic  in  each  of  its  N  variables 
separately  (cf.,  Siegel,  1971,  Vol .  II).  It  follows  that  the  motion 
described  by  an  N-band  solution  of  the  KdV  equation  may  be  thought  of 
as  uniform  translation  of  a  point  along  a  straight-line  path  on  an 
N-dimensional  torus.  Because  any  such  solution  has  only  N  degrees 
of  freedom,  rather  than  infinitely  many  as  (4.1)  suggests  a  priori , 
it  has  a  relatively  short  recurrence  time  tnat  may  be  estimated  from 
a  knowledge  of  { w j ^ j - •)  (4.20b).  This  is  apparently  the 

explanation  of  the  relatively  short  recurrence  time  observed 
numerically  by  Zabusky  and  Kruskal  (1965).  [Recently,  Thyagaraja 
(1979)  proposed  another  explanation  of  the  short  recurrence  times  in 
the  periodic  KdV  problem.  His  method  is  based  on  counting  the 
number  of  Fourier  modes,  rather  than  the  number  of  unstable  bands, 
required  to  describe  a  KdV  solution.  His  method  provides  a  rigorous 
upper  bound  on  the  recurrence  time,  but  it  seems  not  to  give  an 
accurate  estimate  of  the  true  recurrence  time.] 

In  an  abstract  sense,  this  theory  of  the  periodic  KdV  problem 
(or  its  generalization  by  McKean  and  Trubowitz,  1976)  is  complete. 

In  a  practical  sense,  it  is  difficult  to  extract  from  the  theory  the 
numbers  that  one  needs  to  make  comparisons  with  experiments.  [The 
difficulty  here  is  simply  that  hyperelliptic  functions  have  not  yet 
been  reduced  to  a  practical  engineering  tool.]  Here  are  some  of  the 
practical  problems  for  which  satisfactory  answers  are  not  yet 
available. 

(1)  Given  an  N-band  solution  of  (4.1),  what  is  its  recurrence  time? 

(2)  Given  appropriate  initial  data  for  (4.1),  estimate  its 
recurrence  time. 


(  ! 
I 


(3)  Zabusky  (1969)  defined  Tfa  to  be  the  time  of  breakdown  of 
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the  solution  of 


ut  +  uux  =  0  (4.21) 

that  evolves  from  u(x,  0)  =  sin  ( 2ttx/L )  .  He  found  empirically 
that  the  recurrence  time  of  the  corresponding  solution  of 

ut  +  uux  +  52uxxx  =  0  (4.22) 


was 


Tr  =  (0.71/5)  Tb.  (4.23) 

Can  this  formula  be  derived  theoretically? 

(4)  "Soli tons"  have  come  to  mean  special  solutions  of  problems 
solvable  by  1ST  on  the  infinite  interval.  Each  soliton  is 
associated  with  a  discrete  eigenvalue  in  that  formulation. 
However,  the  word  was  coined  by  Zabusky  and  Kruskal  (1965)  to 
describe  phenomena  observed  in  the  periodic  problem.  What 
were  the  "soli tons"  they  observed?  How  does  the  band 
structure  of  a  solution  of  the  periodic  problem  relate  to  the 
number  of  soli  tons  that  one  would  observe  numerically? 
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Part  V 

Deterministic  and  Chaotic  Models 


In  the  course  of  these  lectures,  I  have  tried  to  survey  the 
theory  of  nonlinear  evolution  equations  ("soliton  theories")  that  can 
be  solved  by  1ST.  If  one  were  asked  to  identify  the  single  feature  of 
these  equations  that  identifies  them  as  special,  one  might  choose  the 
predictability  of  their  solutions  as  their  common  identifying  feature. 
This  property  may  be  stated  in  several  ways. 

i)  These  equations  are  deterministic,  not  only  for  some  finite 
time,  but  for  all  time. 

ii)  Their  solutions  are  (neutrally)  stable  with  respect  to 
perturbations  in  the  initial  data. 

iii)  Their  solutions  are  predictable. 

All  of  these  notions  are  closely  related  to  each  other,  and  they 
amount  to  saying  that  for  these  special  equations,  the  initial-value 
problem  makes  sense  in  a  practical,  real-world  way.  Because  of  their 
predictability,  these  completely  integrable  equations  may  be  regarded 
as  prototypes  of  deterministic  problems. 

As  models  of  physical  systems,  they  contrast  sharply  with 
stochastic  models,  which  require  a  certain  amount  of  unpredictability. 
This  unpredictability  may  come  from  a  variety  of  sources;  one 
possibility  is  that  the  dynamics  themselves  may  be  unpredictable.  An 
example  of  this  inherent  unpredictability  may  be  seen  in  a  system  of 
coupled  ordinary  differential  equations  known  as  the  "Lorenz  model": 

i  *  o(y  -  x)  , 

y  +  y  -  -x(z-r)  ,  (5.1) 

z  *  xy  -  bz  , 
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where  (")  =  d(  )/dt  ,  and  (a,  r,  b)  are  parameters.  In  their 
original  context  these  equations  modelled  a  problem  in  convection,  and 
o  represented  the  Prandtl  number,  r  the  Rayleigh  number,  and  b  was 
a  measure  of  a  length  scale.  For  this  discussion,  however,  we  may 
regard  them  simply  as  a  dynamical  system  with  parameters. 

Lorenz  (1963)  chose  certain  values  of  the  parameters 
(a  *  10,  b  =  8/3,  r  =  28),  and  integrated  the  equation  numerically. 
His  results  are  described  in  detail  in  his  original  work,  and  may  be 
summarized  by  saying  that  the  solutions  he  computed  were  unpredictable. 
These  equations  have  become  a  popular  model  of  a  dynamical  system 
exhibiting  chaotic  behavior.  Based  on  a  great  deal  of  numerical 
work  following  that  of  Lorenz,  the  current  consensus  is  that  the 
solutions  of  these  equations  are  ergodic,  although  no  proof  has  been 
found.  [The  original  notion  of  an  ergodic  trajectory,  due  to 
Boltzmann,  was  that  in  the  course  of  time  it  should  wander  everywhere 
in  the  available  phase  space.  Some  technical  revisions  have  been 
required  since  then,  but  the  original  notion  still  is  close  to  being 
correct.  For  a  more  complete  account,  see  Arnold  and  Avez  (1968) 
or  Arnold  (1978).] 

We  may  regard  the  Lorenz  model  with  the  particular  parameters 
chosen  by  Lorenz  as  a  prototype  of  a  "chaotic  model”.  Its  features 
may  be  characterized  in  the  following  ways. 

i)  The  equations  are  deterministic  for  a  finite  time  (i.e.,  given 
finite  initial  data,  a  unique  solution  exists).  However,  the 
solution  becomes  less  and  less  determined  by  the  initial  data 
as  time  increases  without  bound. 

ii)  The  solution  is  unstable  with  respect  to  perturbations  of  the 
initial  data. 

lit)  Given  any  initial  data,  the  solution  is  unpredictable  over  a 
long  time  scale  in  any  practical  sense. 
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Now  we  have  two  kinds  of  dynamical  systems.  The  Lorenz  model  is 
chaotic,  and  one  might  imagine  using  statistical  methods  to  analyze 
the  problem  further.  Soliton  theories  are  deterministic,  and 
statistical  methods  would  only  obscure  matters  there.  Thus,  a 
fundamental  question  is  “How  does  one  know  whether  a  given  model  will 
be  chaotic  or  deterministic?" 

This  same  question  arises  in  the  foundations  of  statistical 
mechanics,  and  of  the  theory  of  turbulence.  Perhaps  the  first  context 
in  which  the  question  arose  was  theology.  In  that  context, 
deterministic  theory  was  called  “predestination",  while  evidence  of 
unpredictable  behavior  was  attributed  to  "free  will".  A  corresponding 
debate,  with  somewhat  different  names,  is  now  in  progress  in 
psychology. 

No  attempt  will  be  made  here  to  classify  all  dynamical  systems  on 
the  basis  of  how  chaotic  their  solutions  are.  We  will  examine  a 
narrower  question,  which  may  be  regarded  as  a  first  step  in 
constructing  such  a  grand  classification  scheme. 

What  determines  whether  a  given  partial  differential 
equation  can  be  solved  by  1ST? 

This  question  has  been  the  motivation  of  recent  work  by  Ablowitz, 
Ramani  and  Segur  (1978,  1980a, b).  The  answer  seems  to  be  related  to 
what  we  may  call  the  "PainlevS  property"  (which  will  be  defined 
shortly) . 

After  this  rather  long  introduction,  we  may  finally  outline  the 
lecture.  First,  we  must  define  the  PainlevS  property  for  ordinary 
differential  equations  (ODEs),  because  everything  else  follows  from 
it.  Next  we  may  show  that  ODEs  with  the  Painlevd  property  are 
related  to  evolution  equations  solvable  by  1ST.  Our  conjecture  about 
characterizing  these  nonlinear  partial  differential  equations  (PDEs) 
then  is  almost  obvious:  they  must  reduce  to  ODEs  of  PainlevS  type. 
What  is  less  obvious  is  how  to  prove  the  conjecture,  but  a  partial 
proof  is  available.  The  notion  that  ODEs  of  P-type  are  closely 
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re  la  ted  to  1ST  and  complete  integrability  may  be  used  in  a  variety  of 
ways.  We  will  look  at  two:  (i)  to  test  whether  a  given  PDE  can  be 
solved  by  1ST;  and  (ii)  to  find  conditions  under  which  the  Lorenz  model 
is  completely  integrate. 


1. 

The  Painlevg  Property 

Consider  first  a  linear  ordinary  differential  equation,  say  of 
second  order: 


■j-y  +  p(z)  *  d(z)w  =  0  .  (5.2) 

az  az 

For  suitable  p(z)  and  q(z)  ,  this  equation  may  be  viewed  in  the 
complex  plane,  and  the  singularities  of  the  solution  of  (5.2)  are 
found  by  examining  p(z)  and  q(z)  (e.g.,  Ince,  1956,  Ch.  15).  In 
particular,  the  general  solution  has  two  constants  of  integration, 

w(z;  A,  B)  =  Awi(z)  +  Bw2(z)  ,  (5.3) 

and  the  location  in  the  complex  plane  of  the  singularities  of  w(z) 
do  not  depend  on  A  or  B  .  The  singularities  of  a  linear  differential 
equation  are  said  to  be  fixed,  because  they  do  not  depend  on  the 
constants  of  integration. 

Nonlinear  differential  equations  lose  this  property.  A  very 
simple  example  of  a  nonlinear  ODE  is 

^  +  Wz  *  0  ,  (5.4) 


its  general  solution  is 


w(z;  z0) 


_J _ 

z-z0  ' 


(5.5) 


Here  z0  is  the  constant  of  Integration,  and  it  also  defines  the 
location  of  the  singularity.  This  singularity  is  movable,  because  its 
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location  depends  on  the  constant  of  integration. 

So  linear  differential  equations  have  only  fixed  singularities, 
nonlinear  equations  can  have  both  fixed  and  movable  singularities. 
About  100  years  ago,  mathematicians  asked  the  following  question: 

Q;  Which  nonlinear  ODEs  admit  no  movable  branch  points 
or  essential  singularities? 


Movable  poles  are  allowed,  as  are  fixed  singularities  of  any  kind.  We 
will  refer  to  this  property  as  the  PainlevS-property,  and  equations 
that  possess  it  will  be  said  to  be  of  Painlev6-type,  or  P-type. 

It  turns  out  that  the  only  first  order  equations  with  the 
Painlevd-property  are  generalized  Riccati  equations 

=  p0(z)  +  pi(z)w  +  p2(z)w2  .  (5.6) 


(A  complete  review  of  the  nineteenth  century  work  in  this  field  may  be 
found  in  Ince,  1956,  Ch.  12-14.) 

Painlevd  and  his  coworkers  were  able  to  answer  the  question 
comprehensively  for  second-order  equations  of  the  form 


d2w 

dz2 


(5.7) 


where  F  is  rational  in  dw/dz  and  w,  and  analytic  in  z.  They 
showed  that  out  of  all  possible  equations  of  the  form  (5.7)  only 
50  canonical  equations  have  the  PainlevS  property  of  no  movable  branch 
points  or  essential  singularities.  Further,  they  showed  that  44  of 
these  equations  can  be  reduced  to  something  already  known,  such  as 
elliptic  functions.  That  left  six  equations  that  defined  new 
transcendental  functions,  called  the  Painlevfi  transcendents.  The 
first  three  of  these  are: 
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d2w 
d z1 


2w3  +  zw  +  a 
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d2w 

dz2 
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There  are  three  more. 

The  question  of  which  equations  have  the  Painlevg-property  is 
appropriate  at  any  order,  but  comprehensive  results  are  available  only 
at  the  first  and  second  order. 


2. 

Relation  to  1ST 

A  relation  between  ODEs  of  P-type  and  the  inverse  scattering 
transform  is  formulated  in  the  following: 

Conjecture  (Ablowitz,  Ramani,  and  Segur,  1978) 

Every  nonlinear  ODE  obtained  by  an  exact  reduction 
of  a  nonlinear  PDE  solvable  by  some  inverse  scattering 
transform  has  the  PainlevS-property. 


Here  are  some  examples.  The  Boussinesq  equation 

utt  s  “xx  +  (f)xx  +  }  uxxxx  (5-8) 

is  a  nonlinear  PDE  solvable  by  1ST  (Zakharov,  1973).  An  exact 
reduction  to  an  ODE  may  be  obtained  by  looking  for  a  traveling  wave 
solution: 


u(x,  t)  =  w(x  -  ct)  =  w(z) 

Then  (5.8)  becomes 


(5.9) 


(1  -c2)w"  + 


*  0  , 


(5.10) 
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which  can  be  integrated  twice.  Depending  on  the  constants  of 
integration,  the  result  after  rescaling  is  either 

w"  +  2w2  +  a  =  0  or  w"  +  2w 2  +  z  -  0  .  (5.11) 


The  first  possibility  defines  an  elliptic  function,  whose  only 
singularities  are  poles.  The  second  possibility  is  the  equation  for 
Pj  .  In  either  case,  the  ODE  has  the  Painlevg  property.  So  the  PDE 
solvable  by  inverse  scattering  reduces  to  an  ODE  of  P-type. 

Another  example  is  the  modified  KdV  equation 

ut  -  6u2ux  +  uxxx  =  0  ,  (5.12) 


which  can  be  solved  by  1ST  (Wadati,  1972).  An  exact  reduction  to  an 
ODE  may  be  obtained  by  looking  for  a  self-similar  solution: 


u(x,  t)  -  (3t)*2,3w(z)  ;  z  =  x/ (3t) 1/3  , 

=*>  w'"  -  6w2w'  -  (zw)'  =  0  . 


(5.13) 


This  can  be  integrated  once 

w"  =  2w3  +  zw  +  a  .  Pjj 

Again,  the  ODE  is  of  P-type. 

The  sine-Gordon  equation 

uxt  *  sinu  (5.14) 

can  be  solved  by  1ST  (Ablowitz,  Kaup,  Newell  and  Segur,  1973).  It  has 
a  self-similar  solution 

u(x,  t)  =  f (z)  ,  z  =  xt  (5.15) 

If  we  set  w(z)  *  exp(if)  ,  then 

w"  *  l(w)2  -  i(w)  ♦  n  Pm 


Again,  the  ODE  Is  of  P-type. 
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By  now  the  pattern  is  evident,  and  I  may  simply  state  that  there 
is  a  nonlinear  evolution  equation  solvable  by  1ST  that  reduces  to  . 
and  another  that  reduces  to  Py  .  The  point  is  not  that  the  evolution 
equation  must  reduce  to  one  of  the  six  Painlevd  transcendents,  which 
are  all  of  second  order,  but  that  it  must  reduce  to  an  ODE  of  P-type. 
We  have  checked  an  enormous  number  of  examples.  In  every  case 
checked,  PDEs  that  can  be  solved  by  1ST  reduce  to  ODEs  of  P-type  and 
PDEs  that  are  not  solvable  by  1ST  (e.g.,  this  may  be  determined  by 

observing  numerically  that  two  solitary  waves  do  not  interact  like 

soli  tons)  reduce  to  ODEs  that  are  not  of  P-type. 

So  there  is  some  kind  of  relation  between  partial  differential 
equations  solvable  by  1ST  and  ordinary  differential  equations  of 
P-type.  This  relation  can  be  used  to  examine  either  the  ODEs  or  the 

PDEs.  To  see  how  it  helps  in  the  study  of  the  ODEs,  consider  the 

mKdV  equation  and  Pjj  .  Recall  that  the  last  step  of  1ST,  the  inverse 
scattering  part,  goes  as  follows.  F(x,  t)  satisfies  a  linear 
partial  differential  equation 

Ft  +  Fxxx  *  0  (5-16> 

subject  to  some  boundary  and  initial  conditions.  Then  K(x,  y;  t) 
satisfies  a  linear  integral  equation  of  the  Gel' fand-Levi tan-Marchenko 
type. 


K(x,y) 


F(x  +  y)  + 


00  oc 

// 


K(x,  z)  F(z+s)  F(s+y)  dz  ds 


y  >  x 


(5.17) 


Once  K  is  known,  then  q(x,  t)  =  K(x,  x;  t)  satisfies  mKdV: 

qt  -  6q2qx  +  qxxx  =  0  .  (5.18) 

In  the  full  1ST  treatment,  F  depends  on  the  initial  data  of 
q(x,  0)  through  the  direct  scattering  problem.  Here  we  simply  start 
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with  F,  and  force  everything  to  be  self-similar: 

5  =  x/(3t) 13  ,  n  =  y/ ( 3t ) 1/3  , 

F(x,  t)  =  (3t)-1'3  F<C)  .  K{x,y;  t)  =  (3t)*‘ 3  <(;,-)  . 

(5.19) 

Then  (5.16)  becomes  a  linear  ODE 

F'"(5)  -  Uf)'  «  0  (5.20) 

and  a  one-parameter  family  of  solutions  is 

F  ( S  +  ", )  •  rAi(Lp)»  (5.21 ) 


where  Ai ( r )  is  the  Airy  function.  The  integral  equation  (5.17) 
becomes: 


K(5.  n)  -  rAi(LTI)+T  f  J  K{Z'  ;)  ^(^T1)  d;de 

n>5  (5.22) 

The  Airy  function  decreases  rapidly  as  its  argument  becomes  large,  so 
the  integral  term  in  (5.22)  is  very  well-behaved.  Therefore,  it  is 
relatively  easy  to  solve  (5.22)  for  n  ^  6  .  On  n  =  S  ,  the  solution 
of  (5.22)  satisfies  the  self-similar  form  of  mKdV,  vi_z . ,  Pjj  : 

K(£,  O  -  2K3te,a  +  SK(E.E)  .  (5.23) 

(Two  different  proofs  of  this  fact  are  given  by  Ablowitz,  Ramani  and 
Segur,  1978,  1980a.)  The  point  here  is  that  (5.22)  is  an  exact 
linearization  of  Pjj  :  every  solution  of  the  linear  integral  equation 
also  solves  Pjj  .  The  general  solution  of  (5.23)  involves  two 
arbitrary  constants;  the  linear  integral  equation  gives  a  one  parameter 


t- 


(r)  family,  which  includes  all  of  the  bounded  real  solutions  of 
(5.23). 

Next  let  us  sketch  a  partial  proof  of  why  this  test  actually 
works.  Consider  a  linear  equation  of  tne  form 


where  F  vanisnes  rapidly  for  large  values  of  the  argument  and  N 
depends  on  F  .  For  example,  in  (2.16),  N(x,  z,y)  =  F(z  +  y).  In 
(5.22),  we  had 


N(x,  2,  y) 


/ 


F ( z  +  s )  F(s  +  y)  ds  . 


(5.25) 


Other  choices  are  also  possible.  We  want  to  show  that  every  solution 
of  a  linear  integral  equation  like  (5.24)  must  have  the  Painlevg 
property.  Then  if  K  also  satisfies  an  QOE,  the  family  of  solutions 
of  the  ODE  obtained  via^  (5.24)  necessarily  has  the  Painlev6  property 
as  well.  So  the  Painlevg  property  is  not  out  of  the  blue,  it  is  a 
consequence  of  the  linear  integral  equation. 

Very  roughly,  the  proof  goes  like  this  (for  details,  see 
Ablowitz,  Ramani  and  Segur,  1980a;  also  McLeod  and  Olver,  1980). 

i)  F  satisfies  a  linear  00E,  and  therefore  has  no  movable 
singularities  at  all . 

ii)  If  F  vanishes  rapidly  enough,  then  the  Fredholm  theory 
of  integral  equations  applies.  It  follows  that  (5.24)  has 
a  unique  solution  in  the  form; 

00 

K(x,y)  =  F(x  +  y)  +  J  F(x  +  z)  ^px»  dz  .  (5.26) 

x 


where  Di  and  Dr  are  entire  functions  of  their  arguments.  Then  the 


i. 


singularities  of  K  can  only  come  from  the  fixed  singularities  of  F  , 
or  the  movable  zeros  of  02  ■  But  02  is  analytic,  so  these  movable 
singularities  must  be  poles. 


3. 

Appl ications 

Here  are  two  examples  of  how  the  conjecture  may  be  used. 

Example: 

In  (1+1)  dimensions,  the  nonlinear  Schr'ddinger  equation  is 

iut  ■  uxx  +  a|u!Ju.  (5.27) 

It  can  be  solved  by  1ST  (Zakharov  and  Shabat,  1972).  A  natural 
general  ization  to  (2  +  1)  dimensions  is 

iut  =  V2  u  +  a | u | 2  u  .  (5.28) 

We  claim  this  equation  cannot  be  solved  by  1ST,  because  (5.28)  has  a 
similarity  solution  in  the  form 

u(x,y,  t)  =  R(/x*  +  y2  ;  >)  exp(ixt)  ,  (5.29) 

and  the  ODE  for  R(r)  is  not  of  P-type.  So  the  nonlinear  Schrodinger 
equation  is  solvable  in  (1  +  1)  dimensions,  but  not  in  (2  +  1) 
dimensions,  or  in  (3  +  1)  dimensions.  On  the  same  grounds,  we  claim 
that  the  equation  for  water  waves  in  deep  water, 

iut  +  uxx  -  uyy  +  o|u|2u  *  0  ,  (5.30) 

cannot  be  solved  by  1ST. 

Example: 

If  the  PainlevS  property  is  as  closely  tied  to  complete 
integrability  as  we  have  claimed,  it  ought  to  identify  values  of  the 
parameters  for  which  the  Lorenz  model  (5.1)  is  completely  integrable. 


-54- 


Thus,  we  may  ask  whether  the  Lorenz  model  is  ever  of  P-type.  The 
answer  is  that  there  are  exactly  four  choices  of  (z,  r,  b)  for  which 
(5.1)  is  of  P-type. 

i)  c  =  0.  In  this  case  the  equations  are  effectively  linear,  and 
therefore  deterministic.  Certainly  the  solutions  exhibit  no 
chaotic  behavior. 

ii)  j  -  1/2  ,  b  =  1  ,  r  -  0  .  The  equations  have  two  exact 
integrals: 

y2  +  z2  =  A2  exp(-2t)  ,  (5.31) 

x2  -  2  =  B  exp(-t)  ,  (5.32) 

after  which  the  third  integration  may  be  obtained  by  quadrature, 
or  the  solution  may  be  expressed  in  terms  of  elliptic  functions. 
The  solutions  may  be  considered  general i zations  of  the  periodic 
orbits  identified  by  Lorenz. 

iii)  o  «  1  ,  b  -  2  ,  r  *  1/9  .  A  first  integral  is 

x2  -  2z  -  C  exp ( -2 1 )  .  (5.33) 

After  an  involved  change  of  variables,  the  resulting  second-order 
equation  becomes  Pjj  .  Again,  the  problem  is  deterministic. 

iv)  a  *  1/3,  b  *  0,  r  arbitrary.  (The  analysis  of  this  case  is 
due  to  A.  Raman i.)  We  may  write  y  *  3x  +  x  from  (5.1a),  and 
replace  (5.1)  with  a  third-order  equation  for  x  .  It  has  a 
first  integral: 

xx  -  x2  +  *  C  exp(-  4/3  t)  .  (5.34) 

With  the  substitution, 

T  =  exp(-t/3)  ,  x ( t )  *  TW(T),  (5.35) 
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(5.34)  reduces  to  Pjjj  with  a  =  6  s  0  .  Thus  the  Lorenz  model 
(5.1)  has  at  least  one  integral,  and  reduces  to  a  classically 
known  equation  of  lower  order,  whenever  the  coefficients 
(c,  r,  b)  are  chosen  so  that  (5.1)  is  of  P-type. 

These  isolated  points  of  parameter  space  are  embedded  in  larger 
regions  in  which  the  equations  have  first  integrals,  although  they 
may  not  be  completely  integrable. 

v)  If  b  =  1  ,  r  =  0,  then  (5.31)  obtains  for  any  o.  The 
existence  of  this  integral  precludes  ergodic  trajectories. 

vi )  If  b  =  2o  ,  then  for  any  (r,  c)  , 

x2  -  2az  *  C  exp(-2ct)  .  (5.36) 

Again,  ergodic  trajectories  are  impossible. 

These  different  possibilities  are  shown  in  a  map  of  parameter 
space  in  Figure  8.  The  point  here  is  that  although  the  Lorenz  model  Fg.8 

may  be  chaotic  for  some  values  of  the  parameters,  it  is  completely 
integrable  for  others,  and  looking  for  the  PainlevG  property  provides 
an  effective  means  of  identifying  points  of  deterministic  behavior. 

Of  course,  this  notion  is  not  restricted  to  the  Lorenz  model. 
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Figure  Captions 

1.  Physical  configuration,  showing  notation  for  (1.9). 

2.  Typical  long-time  solution  of  (2.1),  with  the  wave-train  spreading 
to  the  left.  Modulations  of  the  wave-train  are  determined  by  the 
initial  data. 

3.  Soliton  solution  of  the  KdV  equation  (2.7). 

4.  Typical  long-time  KdV  solution.  The  solitons  are  determined  by 
the  discrete  spectrum  [<(  and  *i  in  this  case],  while  the 
radiation  is  associated  with  the  continuous  spectrum  [b(k)]. 

5.  Measured  packet  of  monochromatic  (frequency  =  1  Hz)  surface  waves 
of  small  amplitude.  The  theoretical  shape  of  the  appropriate 
envelope  soliton  solution  of  (3.16)  is  given  by  the  dashed  line. 
From  Ablowitz  and  Segur  (1979). 

6.  Cnoidal  wave  solution  of  the  KdV  equation  (4.1)  with  k  -  1/2  and 


7.  Floquet  diagram  for  a  particular  u(x)  in  (4.7).  In  this 
example  there  are  four  unstable  bands. 

8.  Map  of  parameter  space  for  Lorenz  model  (5.1),  showing  where 
equations  have  Painlevl  property  (P^  -  P^)  and  where  they  admit 
exact  integrals. 
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